Quantum entanglement of identical particles 
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We consider entanglement in a system of fixed number of identical particles. Since any operation 
should be symmetrized over all the identical particles and there is the precondition that the spatial 
wave functions overlap, the meaning of identical-particle entanglement is fundamentally different 
from that of distinguishable particles. The identical-particle counterpart of the Schmidt basis is 
shown to be the single-particle basis in which the one-particle reduced density matrix is diagonal. 
But it does not play a special role in the issue of entanglement, which depends on the single-particle 
basis chosen. The nonfactorization due to (anti)symmetrization is naturally excluded by using the 
(anti)symmetrized basis or, equivalently, the particle number representation. The natural degrees 
of freedom in quantifying the identical-particle entanglement in a chosen single-particle basis are 
occupation numbers of different single particle basis states. The entanglement between effectively 
distinguishable spins is shown to be a special case of the occupation-number entanglement. 
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How does one characterize entanglement in a fixed 
number of identical particles? Obviously, a correct char- 
acterization must exclude the nonfactorization due to 
(anti)symmetrization. Here, we clarify that it can be 
done by using the (anti) symmetrized basis, which is 
equivalent to the particle number representation. This 
naturally leads to the use of occupation numbers of dif- 
ferent single-particle basis states as the (distinguishable) 
degrees of freedom in quantifying identical-particle en- 
tanglement even when the number of particles is con- 
served. The occupation-numbers of different modes have 
already been used in quantum computing lj. The use 
of modes was made in a previous study of identical- 
particle entanglement, based on formally mapping the 
Fock space to the state space of qubits or harmonic os- 
cillators 0, but it was under the unphysical presump- 
tion of full access to the Fock space. We shall elaborate 
that the concept of entanglement in a system of identical 
particles is fundamentally different from that of distin- 
guishable particles, for which entanglement is invariant 
under local unitary transformations. There is no local 
operation that acts only on one of the identical parti- 
cles. The single-particle basis transformation is made on 
each particle and chooses a different set of particles in 
representing the many-particle system. Thus, the en- 
tanglement property of a system of identical particles 
depends on the single-particle basis used. The parti- 
cle number basis state for a fixed number of particles 
is just the normalized (anti)symmetrized basis in the 
configuration space, i.e., Slater determinants or perma- 
nents. Therefore, the occupation-number entanglement 
in a fixed number of particles is nothing but the situation 
that the state is a superposition of different Slater de- 
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tcrminants or permanents. Another consequence is that 
the two-identical-particle counterpart of the Schmidt de- 
composition, which we call Yang decomposition since the 
corresponding transformation of an antisymmetric ma- 
trix was first obtained by Yang long ago |3|, does not 
play a similar role in characterizing the entanglement. 
On the other hand, we show that like the Schmidt ba- 
sis, the Yang basis is the single-particle basis in which 
the one-particle reduced density matrix is diagonal. It 
is a common wisdom to treat the entanglement between 
spins of identical particles, when they are effectively dis- 
tinguished in terms of another degree of freedom, in the 
way of distinguishable particles. We show that it is in 
fact a special case of the occupation-number entangle- 
ment with a constraint on the accessible subspace of the 
Fock space. 

In terms of the product basis \k±, ■ ■ ■ , kjf) = • -® 

|fcjv), the TV-particle state is 

\i>) = ^2 q(ki,---,k N )\ki ) ---,k N ), (l) 

fel,"-,fcjV 

where the summations are made over k\ , • • • , k^ indepen- 
dently, the coefficients q(k±, • • • , fejv) are (anti)symmetric. 

It is often convenient to use the unnormal- 
ized (anti)symmetrized basis, |fci, • • • , fc/v)^ = 
^2p'(— 1) p |A:i, • • • , fcjv), where P denotes permutations, 
'+' is for bosons while '-' is for fcrmions. Suppose that 
in ki, ■ ■ ■ , fcj\r, there are n a fcj's which arc a, then there 
are only TV!/J|^ n Q ! different permutations. Hence, 
|fci, • • • ,/c N ) (±) = J2'p(^ 1 ) P I\ a n a\ k ir ■ ■ ,k N ), where 
the summation is only over all different permutations. 
The TV-particle state is then 

IV;) = ]T g(k 1 ,---,k rr )\k 1 ,---,k N )^\ (2) 

(fci, 

where (fci, • • • , /cat), disregarding the order of kj_, ■ ■ ■ , fcjv, 
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is a single index. Up to the sign depending on the or- 
der of k\, ■ ■ ■ , in q(ki, • • • , fejy), 9 is equal to g, i.e., 
each set of (anti)symmetrized terms in Eq. corre- 
sponds to one term in Eq. J2J. Equation J2J can be 
rewritten in terms of the normalized (anti)symmetrized 
basis \h,---,k N )^ = y/l/N\Y[ a n a \\k lr --,k N )^\ as 

|V>) = E,(k u -,k N ) h ( k U ■ ■ ■ ,k N )\ki, ■ ■ ■ ,fciv) (s) ■ 

For a fixed number of particles, the normalized 
(anti)symmetrized basis can be rewritten in terms of 
the occupation numbers of different single-particle ba- 
sis states. This is the particle number representation, in 
which 

IV') = ^2 f(ni,---,n 00 )\ni ) ---,n 00 ), (3) 

ni , ■ ",7100 

where rtj is the occupation number of mode j, 
\ni, - ■ ■ ,noo) = (a\) ni • ■ ■ (atc) n °° |0) , the summations are 
subject to the constraint J2 a n a — N, hence in the com- 
plete summationjOf course most of the /'s are zero. 

In Refs. 0, IE IE S 13, an arbitrary iV-particle state, 
in an arbitrary single-particle basis, is inappropriately 
written as £ n ... lJV w ll ... lN a\ i • ••aJjO), where w n ... lN is 
(anti)symmetric, and each subscript of the creation op- 
erators runs over all the modes. One should note that 
the creation or annihilation operators are associated 
with (anti) symmetrized basis. For example, a\a)- 10) = 

±otoj|0) = \li)\lj) = ± \i)\j)), where i + j. 

Therefore a\ ■ ■ ■ a\ N |0) in these papers may be corrected 
to | z i - - - in). On the other hand, if one uses the particle 
number basis states, no (anti)symmetrization needs to 
made Q. 

Single-particle basis transformation for identical parti- 
cles is not the counterpart of the local unitary transfor- 
mation in a system of distinguishable particles. It acts 
on each identical particle in the same way. For distin- 
guishable particles, local unitary transformations do not 
change the entanglement. In contrast, for identical par- 
ticles, the entanglement depends on which single-particle 
basis is chosen. Consequently, unlike the Schmidt decom- 
position of distinguishable particles, the Yang decompo- 
sition does not play a special role in identical-particle 
entanglement (lO|. 

As a simple example, consider a two-particle state 
•WJO) = 75 (l k i>l k 2) ± |k 2 )|ki)), assuming ki ± k 2 . 
In terms of momentum basis, this is only a basis state 
in particle number representation. Written in terms of 
the product basis, the non-factorization is only due to 
(anti)symmctrization. Therefore, there is no entangle- 
ment. However, in terms of the position basis, it becomes 
£ ri ,r, e i ( kl ' ri+k2 r2 )4 i atjO), which is entangled. 

Since a single-particle basis transformation is made on 
every particle, the entanglement property depends on the 
single particle basis chosen. This invalidates the use of 
the von Neumann entropy of the one-particle reduced 
density matrix as a measure of entanglement for two 
identical particles 0- Of course, the von Neumann en- 



tropy of any density matrix characterizes its decomposi- 
tion in its eigenbasis, which is the Yang basis in the case 
of one-particle reduced density matrix for a system of 
two identical particles, as explicitly shown below. The n- 
particle reduced density matrix for a .ZV-particle system is 
(k'ir ■ ■ ,k' n \P {n) \ k i, - ■ ■ , k n) = Tr{a h i t ■ ■ ■ a k > n pa\ n ■ ■ ■ 4j, 
with Trp(") = N{N - 1) • • • (N - n + 1). One can find 

(k[--.k' n \p^\k 1 ---k n ) = w ^ l X 

T,k n+1 -k N (±) ( k 'i ■ --Kkn+i ■ ••fcj V |p|fci • --k n k n+1 ■ ■•A; A r) (±) . 

(4) 

For a two-boson product state aj^aj^ |0) with ki ^ 
k 2 , the one-particle reduced density matrix is given 
by (kilp^lkx) = (k 2 |pW|k 2 ) = 1 and (k^ 1 )^) = 
(k 2 1 p^ 1 ^ |k x ) = 0, hence the one-particle partial entropy 
is log 2, contradicting the previous claim. 

The dependence of entanglement on the single-particle 
basis is consistent with the point of view that individ- 
ual particles are excitation of quantum fields, and that 
each different single-particle basis, in fact, defines a dif- 
ferent set of particles representing the many-body state. 
In fact, in many-body physics, it is a routine to make var- 
ious transformations, which usually changes the nature 
of entanglement • 

With (anti)symmetrization already made on the ba- 
sis, the correlation embedded in the coefficients natu- 
rally gives the information on entanglement. A Slater 
determinant or permanent is just a (anti)symmetrized 
basis state, hence is nonentangled with respect to the 
given single-particle basis. One can regard a superpo- 
sition of Slater determinant or permanent as entangled 
in the given single-particle basis. Any operation, even a 
one-body one, of which the single-particle basis transfor- 
mation is an example, acts on all the particles. A trans- 
formation from a superposition of Slater determinant or 
permanent to a single Slater determinant or permanent, 
in another single-particle basis, must involve operations 
on all particles and actually chooses a different set of 
particles in representing the state. In a sense, there is a 
builtin nonseparability, based on both the symmetriza- 
tion of any operation and the spatial wave function over- 
lap. Consistently, without spatial wave function ove rlap 
or under the condition of the so-called remoteness |12l| . 
the symmctrization does not have any physical effect. 

Hence, the entanglement is between different single- 
particle basis states. Whether a certain single particle 
basis state is entangled with other single particle basis 
states can be decided by whether the former is mixed 
with the latter in the single-particle basis transformation 
which transforms the superposition into a single Slater 
determinant or permanent. This can be seen most clearly 
by using the second quantization. For example, in a two- 
particle state l/Vmfllti^ks + ■ • • + a k m+1 )|0), where kj's 
are different from each other, m > 1, 1^) state is ob- 
viously separated from the others. One can obtain the 
one-particle partial entropy as log 2 + 1/2 log to > log 2. 

Since, the distinguishable label is the set of oc- 
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cupation numbers of different single-particle basis 
states, clearly they can be used to quantify entan- 
glement, in the way of entanglement between dis- 
tinguishable objects. From Eq. J3J), one obtains 
the density matrix as (n' 1} ■ ■ ■ , n'^plrii, ■ ■ ■ , n.^) = 
f*(n'i, • • • , n' 00 )/(ni, • • • , rtoo), from which one can ob- 
tain the reduced density matrices of occupation 
numbers. For example, the reduced density ma- 
trix of mode 1 is defined as (n' 1 |pi(l)|rii) — 
E„ 2 ,., nco (n'x, «2> • • • ,noo\p\n>x,ri2 ■ ■ ■ , n^). Similarly, 
the reduced density matrix of the set of modes 1, • • • , / 
is 

(n' 1 ,---,n' l \pi(l---l)\ni,---,ni) = 

T,n l+1 ■■■ n ' ' •j^j^+Ij^ooIpK, • ■ ■ , 71;, 7lj+i, 

(5) 

the nonvanishing elements of which satisfy X)i=i n i = 
y~li—i n% as constrained by the particle number conserva- 
tion. From these Fock-space reduced density matrices, 
one can, for example, calculate bipartite entanglement 
between the occupation numbers of I modes and the oc- 
cupation numbers of the other modes. 

It is important to note that the use of occupation num- 
bers as the degrees of freedom in characterizing entangle- 
ment is valid even when the particle number is conserved. 
This physical constraint, as well as the constraints that 
for fermions rij is either or 1 and that the number of 
the relevant modes [D| may be finite are all automati- 
cally satisfied by the set of nonzero f's. Hence, this ap- 
proach is a natural one within the standard second quan- 
tization formalism, compatible with the representations 
of the observables in terms of creation or annihilation 
operators, which can be viewed as coordinated transfor- 
mations of occupation numbers of a set of modes. The 
second-quantized representation of an n-body operator 
O is J2 a \ ' ' ' a l (*i ' ' ^'JOlh ■ ■ -i n )ai n ■■■ai 1 - One can 
observe that, for example, there is no operation which 
only changes the occupation number of one mode. One 
may consider "second-quantized computation" . 

In principle, one can define entanglement with respect 
to any reference state of the system. In this case, the 
occupation number of each mode in defining the relative 
entanglement is the difference with that in the reference 
state, as conveniently seen by considering the action of 
creation operators. There are two reference states that 
are of particular interest. One is the empty state, as we 
have implicitly considered up to now. Another one is 
the ground state of the system, which is suitable when 
all physical processes are in a same bulk of material. In 
discussing the entanglement in a ground state, it is with 
respect to the empty state. The considerations can even 
be extended to relativistic quantum field theory, where 
the ground state is the vacuum. 

An important situation is that the single-particle basis 
includes both spin and orbit (momentum or position). 
One can denote the total index as (o, s), where o sub- 
stitutes for momentum k or position r. A special case 
is half-filling, i.e., each orbit is constrained to be occu- 



pied by only one particle, i.e. ^ s n o s = 1 for each rele- 
vant value of o. Then with the orbit modes as the labels 
with which the particles are effectively distinguished, the 
entanglement can be viewed as the spin entanglement 
among the particles in different orbit modes. Under the 
constraint of half filling, the many-particle state must 
be a sum of products of Y[ s \ n o,s), under the constraint 
^2 s n o,s = 1, for relevant orbits. In the many-particle 
state, one simply rewrites JT S |ra ,«) as I So), where So is 
unambiguously the s corresponding to n 0jS — 1. This 
rigorously justifies the common wisdom that although it 
is meaningless to identify which particle is in which or- 
bit, it is meaningful to say that the particle in a certain 
orbit is spin entangled with the particle in another or- 
bit. Hence, the entanglement between Heisenberg spins, 
which appears as entanglement between distinguishable 
objects, is in fact a special case of occupation-number 
entanglement. 

Entanglement between Heisenberg spins is the basis 
of the quantum computing scheme based on electrons 
in double quantum dots [l4[- When the electrons are 
separated in the two dots, because there are only one- 
dot potentials, while the Coulomb interaction is negli- 
gible, the condition of remoteness ^2 is satisfied. One 
can verify that an antisymmetrization between electrons 
in different dots has no physical effects. On the other 
hand, when they are close, and the interaction is appre- 
ciable, the antisymmetrization has physical effects, and 
the entanglement can be characterized by using the full 
formalism of occupation- number entanglement. During 
the interaction period, they access the full Hilbert space, 
which includes the state in which the two electrons, with 
opposite spins, locate in a same dot, i.e. |l)i,f)|li,|)> 
where i — 1, 2 represents the dots. There are 4!/2!2! = 6 
two-particle antisymmetrized basis states, or occupation- 
number basis states. Nevertheless, as far as the ini- 
tial and final states are in single occupancy, and the 
Heisenberg model is a valid description, the interaction 
period can be viewed as an intermediate process deter- 
mining the effective spin coupling, while the leakage into 
the full Hilbert space of two identical particles during 
a two-particle gate operation does not cause any prob- 
lem. In terms of the occupation-number states, the spin 
state of each electron in each dot is | |)i = |l)t,t |0)i,l, 
| J.), = 10)^1-11)^1. Because a spin qubit is in fact an 
occupation-number state, the loss of identification after 
separating from the double occupation, as concerned in 
Ref. 0|, does not matter. Note that the intermediate 
state with double occupancy is necessary for the electrons 
to interact in order to undergo a two-qubit operation. 

Finally, we come to the question that what is special 
about the Yang basis. For two distinguishable parti- 
cles, the Schmidt basis is clearly the one in which the 
reduced density matrix of each particle is diagonal: For 
J2i c iV)a\i)b of distinguishable particles a and b, the ele- 
ments of the reduced density matrix of either a or & are 
given by (i\p a (b)\j) = l c i| 2 ^j- I n t ne following, we show 
that the Yang basis is the basis in which the one-particle 
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reduced density matrix is diagonal. 

In their Yang basis, a two-fermion state is like \ipf) — 
ci(|l)|2)-|2)|l)) + C2(|3)|4)-|3)|4)) + --- where we use |») 
to denote different single-particle basis states. If k 2 = k2, 
then k' x — k\ is necessary for any of (k[k' 2 \p\kik2), 
(k' 1 k' 2 \p\k2ki), (k' 2 ki\p\kik2) and {k' 2 k\ \p\k2k1) to be 
nonvanishing. Therefore using Eq. one finds 

(fci|pW|fci) = S k[kl J2 k2 C-><fcifc2|plfci*?a> c ->. Hence, p« 
is diagonal. 

In their Yang basis, a two-boson state is like ji/'b) = 
di|l)|l) + d 2 |2)|2) + • • Then one finds (fc;^|p|fcifc 2 ) = 
(k'^p^h) = (fc^ilplfcife) = (k'^p^h) = 
Sk' 1 k' 2 Sk 1 k 2 (kiki\p\k2k 2 ) ■ Consequently, using Eq. (QJ, one 

finds (fellp^lM - Hence, p« 

is diagonal. 

Let us summarize. If one uses the prod- 
uct basis, the coefficients mix the information on 
(anti)symmetrization and that on entanglement. If, in- 
stead, the (anti)symmetrization is made on the basis, 
then the coefficients unambiguously give the informa- 
tion on entanglement, with respect to the given single- 
particle basis. (Anti)symmetrized basis is equivalent to 
particle number representation, and the occupation num- 
bers of different modes are distinguishable degrees of free- 
dom which can be used in quantifying the entanglement 



even when particle number is conserved. Entanglement 
of identical particles is a property dependent on which 
single-particle basis is chosen, as any operation should 
act on each identical particle in the same way. Indeed, 
individual particles are excitations of a quantum field, 
and the single-particle basis defines which set of particles 
are used in representing the many-particle state. The 
many-particle state is entangled in the corresponding sin- 
gle particle basis when it is not a single Slater determi- 
nant or permanent. The entanglement is between dif- 
ferent single-particle basis states in the given basis. We 
also show that the entanglement between effectively dis- 
tinguishable spins of identical particles is a special case 
of the occupation- number entanglement. We have dis- 
cussed its use in quantum computing. The (necessary) 
leakage into the larger Hilbert space during the inter- 
mediate two-particle process is harmless. Finally it is 
shown that the two-identical-particle counterpart of the 
Schmidt basis is the basis in which the one-particle re- 
duced density matrix is diagonal. In addition to quantum 
computing implementations involving identical particles, 
the result here is also useful for many-body physics |l6j | . 

This publication is an output from project activity 
funded by The Cambridge-MIT Institute Limited. 
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